Classical dispersion-cancellation interferometry 
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Even-order dispersion cancellation, an effect previously identified with frequency-entangled pho- 
tons, is demonstrated experimentally for the first time with a linear, classical interferometer. A 
combination of a broad bandwidth laser and a high resolution spectrometer was used to measure 
the intensity correlations between anti-correlated optical frequencies. Only 14% broadening of the 
correlation signal is observed when significant material dispersion, enough to broaden the regular 
interferogram by 4250%, is introduced into one arm of the interferometer. 



Interferometry is an indispensable tool for precision mea- 
surements. Low-coherence, or white-light, interferometry 
is used for precise measurements of material properties, 
such as optical path length and dispersion. Optical co- 
herence tomography (OCT), a technique for non-invasive 
medical imaging, is based on low-coherence interferom- 
etry [TJ [2]. Both white light interferometry and OCT 
use broad bandwidth light sources to achieve micrometer 
scale image resolution [3 . Although large spectral band- 
width is essential for obtaining high resolution, it also 
increases dispersive broadening of the interferogram. 

Quantum metrology uses quantum mechanical fea- 
tures, such as entanglement and squeezed light, to im- 
prove the sensitivity of measurement devices [4j [5]. 
A two-photon quantum interferometer [6], based on 
frequency-entangled photon pairs, has been demon- 
strated to be insensitive to all even orders of dispersion 
(dEL also see [9j). This effect, known as quantum dis- 
persion cancellation, was proposed as the basis for quan- 
tum-optical coherence tomography [10] and a proof-of- 
principle was demonstrated experimentally [11]. A very 
recent theoretical model [12] claims that interferomet- 
ric dispersion cancellation does not require the use of 
individual pairs of entangled photons. The scheme is in- 
stead based on a nonlinear optical interferometer that 
employs broad-band phase conjugation between two re- 
flections from the same sample. Experimental implemen- 
tation of this technique would be extremely difficult re- 
quiring both development of novel optical sources and a 
suitable method of phase conjugation. Other approaches 
use numerical methods to compensate dispersion in data 
and images obtained with low-coherence interferometry 

or oct nam[i3[i![i7]. 

Dispersion cancellation is staight forward in quantum 
interferometry, but the methods proposed so far in clas- 
sical interferometry are not. Can we use the intuition 
derived from quantum technologies to achieve dispersion 
cancellation in a simpler way in a classical interferome- 
ter p~8j [19] [20]? In this work, we show that dispersion 
cancellation can be achieved using only a classical light 
source, linear optics, and frequency-correlated detection. 



We review quantum dispersion cancellation and use its 
essentials to design an analogous classical system. 

Consider the nonclassical two-photon interferometer 
shown in Fig. la) [6]. The upper path is of length, Li, 
and the lower path is of length, L2=L\ + A. A non- 
linear crystal, pumped by a narrow bandwidth laser of 
frequency 2o;o, generates photon pairs with central fre- 
quency ojq, via parametric down- conversion into the up- 
per and lower paths of the interferometer in the state, 
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The subscripts 1 and 2 are mode labels, and A(ujq-\-5uj) is 
the amplitude for a pair of photons of frequencies ujo + Suj 
in mode 1 and c^o — Su in mode 2. The sum of the fre- 
quencies in each term of the superposition is fixed by 
energy conservation - this is a frequency-entangled state 
with perfect frequency anti-correlation. The photons are 
interferometrically combined at the 50/50 beamsplitter 
followed by two single-photon counting detectors. The 
signal of interest is the number of coincident photon de- 
tection events as a function of the optical delay, A. 

Insertion of a dispersive, lossless medium of length, 
L, in the upper path of the interferometer results 
in a frequency-dependent phase shift, (/>m(^)=&m(^)£, 
where &m(^) is the wavevector at frequency, oj, in the 
material. We series expand &m(^) about ujq- 
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The first derivative is the inverse of the group velocity, 
v g , at ujo and leads to a shift in the centre of the interfer- 
ence pattern. The second derivative is the leading-order 
dispersive term, which causes loss of both spatial reso- 
lution and contrast in low-coherence interferometry by 
broadening the width and reducing the visibility of the 
interference pattern. 

Following Ref.[7 , we make the assumption A(ujq + 
Suj)=A(ujo — Suj), i.e., the amplitude is symmetric about 
the central frequency c^o- We find the expected coinci- 
dence rate, as a function of A is: 



2 



C(A) 



/ 



dduj \A{ujq + Sco)\ 2 
2Sco(L + A) 



1 — cos 



2 Aa, 



(3) 



The expression is in agreement with Ref. [7 j, but uses 
slightly different notation. Notice that the second deriva- 
tive does not appear - this is the dispersion cancellation. 
In fact, all higher-order even derivatives are cancelled. 
Maximum destructive interference occurs when the ar- 
gument of the cosine term is zero for every frequency. 
This happens when the extra group delay imposed by 
the material is balanced by extra optical delay in the 
other arm of the interferometer. We will refer back to 
this expression when describing our classical system. 

Now consider the Mach-Zehnder interferometer in 
Fig. lb). The dimensions and mode labels of this inter- 
ferometer are identical to that described in Fig. la, as is 
the dimension of the dispersive material; both beamsplit- 
ters are 50/50. The intensity spectrum of the input light 
is I (w). The intensities registered by the spectrometers 
for a delay position, A, and frequency, a;, in the outputs 
labeled a and b are, 



I a ((j,A) = /(cj)cos 2 
I b (u,A) = I (uj) sin 2 
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Each of these intensities is affected by all orders of dis- 
persion in the series expansion of <j>M (^)- 

Quantum dispersion cancellation is a result of 
frequency-entanglement in fourth-order, i.e., coincidence, 
detection. Our approach seeks to mimic this effect as 
closely as classical physics allows. We use frequency cor- 
relations, the classical analogue to entanglement, and 
measure a fourth-order signal, achieved by multiplying 
pairs of intensity measurements. Specifically, we mea- 
sure the signal, «S, 
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The integrand of this function is the product of two 
intensities with an energy sum of 2uoq. We use Eqs. 
[4] & [5] and assume, as we did in the quantum case, 
that the input spectrum is symmetric about c^o, i.e., 
I(uoq + 5) = I(uoq — 5), to obtain, 



5(A) = \jdfiw [/(w 



■ 5co)Y 



(7) 



COS 



COS 



2a; (L+A) 



2Lk 



This is the signal of interest from our classical system 
and can be directly compared with the quantum signal in 
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FIG. 1: Dispersion cancellation interferometry in a a) Two- 
photon interferometer using frequency-entangled photon pairs 
[3 [8] and b) a white-light Mach-Zehnder interferometer with 
frequency-correlated detection. These interferometers and 
their expected output signals are described in the text. Co- 
inc. is coincidence detection; S describes a function of the 
output from the spectrometers, c) Experimental realization. 
A broadband laser is the source for a fibre-based two-path 
(Michelson) interferometer. The setup uses a 50/50 beam- 
splitter (FC 50:50), polarization controllers (PC), collimat- 
ing lenses (CL), neutral-density filters (F), two BK7 prisms 
for dispersion control (DC), a translation stage (TS), mirrors 
(M), and a spectrometer. The spectrometer contains a CL, a 
diffraction grating (DG), and focusing lens (FL). 

Eq|3} The argument in the first cosine term is identical 
to the quantum expression and describes a dispersion- 
cancelled interference dip. The second cosine term does 
not appear in the quantum case. Notice though, that its 
argument has only weak dependence on the frequency dif- 
ference 5uo (the integration variable) through the disper- 
sion term. It describes a rapidly oscillating component, 
with wavelength Ao = ttc/uq, with a slowly decaying en- 
velope. The separation of length-scales between these 
terms allows removal of the unwanted fast oscillation in 
the final data with for example, a low-pass Fourier filter. 
The other feature of interest in our classical expression is 
that the first term is multiplied by 1/2. This imposes the 
well-known classical limit of 50% on the destructive in- 
terference visibility in two-photon interferometers [6j [21] . 
The signal S is the classical analogue to the Hong-Ou- 
Mandel dip [6] and contains the same resistance to dis- 
persion as its quantum counterpart. 

The experimental setup is shown in Fig. 2. A compact, 
fiber-pigtailed, femtosecond laser (Femtolasers Inc., cen- 
tre wavelength 792nm, bandwidth FWHM 154nm, aver- 
age power 60mW) was coupled to a fiber-based Michel- 
son interferometer. Broad bandwidth optical and fiber 
optic components were chosen to support propagation 
of the entire laser bandwidth with minimal spectral and 
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power losses. A pair of BK7 prisms mounted on minia- 
ture translation stages in the reference arm of the system 
were used to precisely compensate material dispersion 
mismatch between the two arms of the interferometer. 
The focusing lens and the mirror in the reference arm 
of the system were mounted on a computer-controlled 
translation stage for variable optical delay. The interfer- 
ence pattern generated by light reflected from the sample 
and reference mirrors was detected with a high-resolution 
(0.09nm) and high-speed (20 kHz readout rate) spec- 
trometer and recorded by a computer. The spectrom- 
eter utilized a 4096 pixel linear- array CCD camera and 
it was calibrated for the spectral range 607nm to 1012nm. 
To demonstrate dispersion cancellation with the clas- 
sical interferometer, measurements were made both in 
a dispersion-balanced system and when flat, uncoated, 
BK7 optical windows of thickness 4.690 ± 0.005mm, 
5.940±0.005mm, and 6.170±0.005mm (and several possi- 
ble combinations) were introduced into the sample arm. 
For each measurement, the reference mirror was trans- 
lated in steps of 0.1/im and the spectral interference 
fringes were acquired with a readout time of 60/is per 
step - at least 4 orders of magnitude shorter as compared 
with typical measurement times in entangled photon ex- 
periments. 

The calculation of the signal function S was performed 
in the following way. One spectrometer reading was taken 
for each motor position to provide us with i a (A, A). The 
wavelength scale was converted to frequency and non- 
linear interpolation was used to extract intensities at 
evenly spaced intervals. We obtained I(oS) by measur- 
ing the intensity from the sample and reference arm sep- 
arately and doubling their sum. Energy conservation, 
I((j)=I a ((j, A)+Ib(u; : A), was applied to extract h((d, A) 
without the necessity for a second spectrometer. To sat- 
isfy the assumption in our theory that I(uo)=I{2uoq — a;), 
I (w) and I a ((d) was multiplied by a mirror version of I ((d) 
with respect to the centre frequency c^o- The integral S 
was approximated by a discrete sum over 4096 equally 
spaced energies. 

The total intensity registered by the spectrometer was 
obtained by adding the intensity measured at each pixel 
at a fixed translation stage position. This signal is equiv- 
alent to a signal that one would have been measured by a 
square-law detector, such as a photo-diode. Two exam- 
ples of the total intensity measured as a function of the 
translation stage position are shown in Fig. 2 for the cases 
where no glass (a) and a 16.800 db 0.009mm thick BK7 
glass window (b) were inserted into the sample arm of the 
interferometer. As a result of the material dispersion, the 
intensity interference pattern is dramatically broadened, 
from (2.04 db 0.03)/im to (88.6 db 0.9)/im, and the fringe 
visibility is reduced [22], from 78% to 14%. This figure 
clearly shows the detrimental effect that dispersion has 
on interference. 

The corresponding correlation signal function, «S, for 
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FIG. 2: Experimental Data, a) & b) Total intensity, as mea- 
sured by summing the intensities measured at each frequency 
by the spectrometer, versus motor position with and two 
passes through (16.800 =b 0.009)mm of BK7 glass in the sam- 
ple arm of the interferometer, respectively, c) & d) S versus 
motor position with no BK7 and 16.8mm of BK7 in the sam- 
ple arm. In e) & f), the data from c) & d) have been subject 
to a Fourier low-pass filter to remove rapidly oscillating terms. 
The solid curves are Gaussian fits. These data show that S 
broadens by only about 14% by addition of the glass while 
the standard intensity interference pattern is broadened by 
4250%. 

the two cases of no material dispersion and 16.800mm 
BK7 glass are shown in Figs. 2c) and 2d), respectively. 
Each of these signals has a sharp dip in addition to a 
rapidly oscillating component that corresponds to the fi- 
nal cosine term in Eqj7[ Note that the magnitude of the 
fast oscillating signal is reduced when a large amount of 
dispersion is present in the interferometer. A similar ef- 
fect was observed when we simulated the measurements 
with a computer model. The data from Fig. 2c) and 2d) 
was filtered with a low-pass Fourier filter to remove the 
fast oscillating term and the the filtered data is presented 
in Figs, e) and f). These dips were fitted to a Gaussian 
function to extract their centres and FWHM. While the 
intensity interference pattern is broadened by 4250% of 
its original size due to material dispersion, the correlation 
signal S is broadened only by 14%. The visibilities of the 
correlation signal dip [22] is reduced from (40.8 ±0.14)% 
and (30.0 ± 0.3)% for Figs, e) and f), respectively (re- 
call that the theoretical maximum visibility is 50%). The 
deviations from the Gaussian shape of the fitting func- 



4 




BK7 Glass - double thickness (mm) BK7 Glass - double thickness (mm) 



FIG. 3: Performance of classical dispersion cancellation, a) 
The width of the interference patterns (total intensity, open 
circle; S, solid circle), as measured by translation stage dis- 
placement, versus the thickness of the glass traversed by the 
beam, i.e., twice the glass thickness. The inset expands the 
y-axis to show the almost constant width of S over the whole 
range of glass thicknesses, b) The shift in the centre of the 
interference pattern versus the thickness of the glass. As dis- 
cussed in the text, these data show that the interference pat- 
tern is displaced by the group delay. 

tion are due partially to the non-Gaussian spectrum of 
the laser as well as the present of higher-order material 
dispersion. 

Fig. 3a) is a plot of the total intensity FWHM (open 
circles) and the correlation signal S FWHM (solid circles) 
as functions of twice the physical thickness of the BK7 
optical flats (we use double the thickness because the 
Michelson interferometer uses a reflection geometry). 

The relative shift in the correlation signal dip centre 
as a function of twice the physical thickness of the glass 
in presented in Fig 3b). We estimated the statistical er- 
ror in the centre of the dip to be about 1/mi based on 
the standard deviation of 5 consecutive measurements. 
This is about a factor of 50-100 times larger than the 
fitting error and could be improved by using a higher 
precision translation stage. However, the most signifi- 
cant measurement error is associated with the widths of 
the BK7 flats. We expect the shift in the centre of the 
dip to be determined by the group velocity by the rela- 
tion, (c a i r /v g — 1)L, where c a i r is the speed of light in 
air and v g is the group velocity, in our case at a wave- 
length of 800nm. From Fig 3b), we evaluate the slope 
0.2633. The accepted value, 0.2631, was obtained from 
the Sellmeier equation [23 . Our errors are dominated by 
the uncertainty in the material thickness, which is about 
0.1%, and at this level the slope from the data and the 
theory agree. The centre of the correlation signal S is 
determined by the group velocity. 

We have theoretically derived and experimentally 
demonstrated a method for cancelling even-order disper- 
sion in classical low-coherence interferometry. Dispersion 
cancellation is not a uniquely quantum effect, since it can 
also be observed in completely classical systems. How- 
ever, the interference visibility in our classical analogue 
is only half that achievable in quantum interferometers 
[6] |5]. Two seemingly contradictory constraints are es- 
sential in both the quantum and classical techniques: a 



wide bandwidth of frequencies provides good time resolu- 
tion, whereas narrow frequency correlations reduce sensi- 
tivity to dispersion. Our approach dramatically reduces 
experimental barriers for dispersion cancellation in low- 
coherence interferometry and optical coherence tomogra- 
phy. 
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